clc;
clear;
A=sym([16,2,3,13; ...
        5,11,10,8; ...
        9,7,6,12; ...
        4,14,15,1]);
E=sym(eye(4))
E([1,4],:) = E([4,1],:);
E
A1=E*A,
A2=A*E